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Abstract 

We have derived masses and ages for 79 early-type galaxies (ETGs) in different environments at 
, z ~ 1.3 in the Lynx supercluster and in the GOODS/CDF-S field using multiwavelength (0.6-4.5 /xm; 

KPNO, Palomar, Keck, HST, Spitzer) datasets. At this rcdshift the contribution of the TP-AGB 
phase is important for ETGs, and the mass and age estimates depend on the choice of the stellar 
population model used in the spectral energy distribution fits. We describe in detail the differences 
among model predictions for a large range of galaxy ages, sho'wing the dependence of these differences 
on age. Current models still yield large uncertainties. While recent models from Maraston and 
Chariot & Bruzual offer better modeling of the TP-AGB phase with respect to less recent Bruzual 
& Chariot models, their predictions do not often match. The modeling of this TP-AGB phase has 
a significant impact on the derived parameters for galaxies observed at high- rcdshift. Some of our 
results do not depend on the choice of the model: for all models, the most massive galaxies are the 
oldest ones, independent of the environment. When using Maraston and Chariot & Bruzual models, 
' the mass distribution is similar in the clusters and in the groups, "whereas in our field sample there is 

QhI a deficit of massive (M > IO^^Mq) ETGs. According to those last models, ETGs belonging to the 

cluster environment host on average older stars "with respect to group and field populations. This 
difference is less significant than the age difference in galaxies of different masses. 
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1. INTRODUCTION 

Superclustcrs of galaxies arc the largest structures ob- 
served in the Universe. Their dimensions can range be- 
tween 10 and '^100 Mpc and are composed of two or 
more galaxy clusters and surrounding groups. These 
structures span a large range in galaxy projected number 
density and permit us to study galaxies seen at the same 
epoch but in very different environments. The study of 
superclustcrs at high redshifts gives us deep insight into 
the role of environment in clusters and groups in the very 
early stages of cluster assembly. Deep multiwavelength 
surveys focusing on the study of superclustcr s at z < 1 
(e.g. CL1604 and the ORELSE program; IGal et al.l 
120081 ILubin et ahl |2009|) have shown in detail the vari- 
ations in star formation rates and galaxy populations as 
a function of environments. It has only been in recent 
years t hat superclustcrs have been discovered a t redshift 
z > 1 (|Nakata et al.l 120051 iTanaka et aIll2009D . In this 
paper we study the early-type galaxy (ETC) populations 
in one of those superstructures, the Lynx superclustcr. 

Lynx is a high redshift (z ~ 1.26) superclustcr, com- 
posed of two clusters (RX J0849-I-4452, hereafter Lynx E 
and RX J0 848-I-4453, hereafter Lynx W) and surround- 
ing groups (jNakata et al.|[2005l : Mei et al., in preparation 
- here after Mil). Th e Lynx W cluster was first identi- 
fied bv [Stanford et all (|1997[ ) as an overdensity in a near- 
infrared (NIR) imaging survey, and then spectroscopi- 
cally confirmed at z 



(Group 2), Zs 



1.264 ±0.003 (Group 3) (Mil) 



spec 



1.2 73; it was later detec ted in 
a deep Chandra observation (iStanford et al.ll200lD. The 
Lynx E cluster was initially found bv lRosati 
as a faint extended X-ray source in the ROSAT Deep 
Cluster Surveys (RCDS) and then spectroscopically con- 
firmed at Zspec = 1.261. The two clusters present differ- 
ent structures. Lynx E shows a compact g alaxy distribu- 
tion, with a central b right galaxy merger (jYamada et al.l 
l2002tlMei eralll2006[) . while Lynx W appears to be at an 
earlier stage of assembly, with a lack of a clear cD, and 
galaxies distributed less concentrated. Their X-ray emis- 
sion gives luminosities of L^^' = (2.8 ±0.2) x lO^"*^ erg.s"^ 



^ for Lynx E and Lynx 
rd et al.l[200l 



and L^' = (1.0 ±0.7) x 10^"^ erg.s~^ tor Lyi 

W respectivelv (iRosati et al.lll999l : iStanfo: 

lEttori et al.ll200'fl ). Their spatial distribution, more com- 
pact for Lynx E and more elongated for Lynx W, also 
indicates that L ynx E is likely mo re dynamically evolved. 
Estimates from iJee et ahl (|2006D weak lensing analysis 
indicate velocity dispersions of ct = 740^^34 km.s~^ and 
(T = 762li[3| km.s~^ for Lynx E and Lynx W, respec- 
tively. Those values are consistent with the spectro- 
sc opic measureme nts of c = 720 ± 140 km.s^^ for Lynx 
E (IJee et al.ll2006D and cr = 650 ± 170 km.s;^ for Lynx 
W (IStanford eralll200l (see also lMei et aDl200l Table 
1). 

Around the two clusters, INakata et al.l (j2005[) iden- 
tified seven group candidates, using photometric 
redshifts derived from optical irnaging with the Sub- 
aru /SupnmeC am (|Mivazaki et al.l [20021 ). indicating the 
presence of a superclustcr extending over an angular 
distance of ~ 20', which corresponds to a luminosity 
distance of ~10 Mpc at the clusters' redshift. Optical 
spectroscopy from the Subaru^ Keck and Gemini tele- 
scopes have confirmed three of those group candidates 
as being part of the superclustcr, with redshifts of 
Zspec = 1-266 ± 0.005 (Group 1), z^pee = 1.262 ± 0.005 



Superclustcrs host large populations of ETGs that are 
mostly found in high density regions. Even if ETGs are 
relatively simple galaxies, large bulges dominated by dis- 
persion velocities, their formation and evolution is still 
not well understood. In the local Universe they define 
a ver y tight red sequ ence in the color- n iagnit ude space 
(e.g.. lBernardi et al.ll200l iBaldrv et al.ll2004[) . that be- 
gins to dissip ate only at the highe st redshifts of known 
clusters (e.g.. iKodama et al.ll200'7[) . This suggests that 
galaxy star formation was quenched in the past and bluer 
galaxies migrate to the red sequence. Multi- wavelength 
studies of ETGs on the red sequence in different envi- 
ronments give us constrains o n the galaxy ages an d their 
star formation history (Mil; IRettura et al.ll20"Tol ). 

In this paper, we present a deep, panoramic multi- 
wavelength survey of the Lynx supercluster, ranging from 
the rest-frame ultraviolet to the infrared. We will quan- 
tify the environmental dependences of the ETC ages 
and masses, by comparing cluster and group galaxies 
to a field sample at the same redshift selected from 
the GOODS (Great Obse rvatories Origins Deep Sur- 
vey, iGiavalisco et al.l 120041 ) obse rvations of the Chan- 
dra Deep Field South (CDF-SI (IGiayalisco et al.l l200l 
iNonino et alll2009l: iRetzlaff et al.lboiOl: Dickinson et al, 
in preparation) . 

In §2, we briefly present the observations, the data re- 
duction and the selection of the galaxy sample on which 
our study relies. We describe the photometry in §3 and 
the spectral energy distribution (SED) fitting method we 
use to derive stellar masses and ages in §4. We study in 
§5 the systematics in the SED fitting, especially the in- 
fluence of the chosen model. We then present our results 
and discuss them in §6 and §7. 

In this paper, we adopt a standard cosmology with 
Ho = 70 km.s-i.Mpc-\ rj„ = 0.30 and Qa = 0.70. 
Unless otherwise specifled, all magnitudes are in the AB 
system and have been corrected for G alactic extinction 
using the maps of ISchleeel et"all (119981) . 

2. OBSERVATIONS 

We have obtained images of the Lynx supercluster 
from the optical to the far-infrared (0.6-4.5 /um), in seven 
bandpasses : R, HST/ ACS F775W and F850LP - here- 
after 1775 a-nd zg5o, J, Kg, Spitzer /IRAC chl and ch2 
- hereafter [3.6/im] and [4.5/im]. Information about the 
data are summarized in Tab [1] and Tab 121 

HST Advanced Camera for Surveys (ACS) observa- 
tions of the Lynx superclustcrs have been carried out as 
part of the ACS Intermediate Redshift Cluster Survey 
(Guaranteed Time Observation, or GTO, program 9 919; 
PI: H. Ford; iFord et al.|[200l IPostman et aLl [20051) for 
the two clusters and from an HST /GO program 10574 
(PI: S. Mei) for the groups. The HST ^775 and Zsso 
cluster (group) imaging was carried out with the ACS 
Wide Field Camera in March and April 2004 (December 
2005 to February 2006) for a total exposure time of 7290 
s and 12220 s (6300 s and 10500 s), respectively The 
ACS WFC resolution (pixel scale) is 0.05" and its field 
of view is 210 " x 240 ". The images were pro cessed 
with the APSIS pipeline (jBlakeslee et al.l I2003D . with 
a LanczosS interpolation kernel. We adopted AB 
photometric zeropoints of 25.678 and 24.867 mag re- 



EARLY-TYPE GALAXIES AT Z - 1.3 



3 



spcctivcly in 4775 and zgso from the HST/ ACS wcbsitcFI. 

The R band images come from two different telescopes. 
The clusters have been observed with the Keck telescope, 
and a wider area including Group 1 and Group 2 has been 
observed with the Palomar telescope. Group 3 has not 
been covered by our R band imaging. 

The Palomar R band imaging (P.L: D. Stern) was ob- 
tained in Novemb er 1999 with the COSMIC instrument 
(|Kells et al.lli99ah in 18 exposures of 900 s each for a to- 
tal exposure time of 16200 s. COSMIC has a resolution 
of 0.2468 "/pixel and a field of view of 9.7' x 9.7'. The 
images were reduced using standard procedures: the im- 
ages were bias corrected, trimmed and flat fielded using 
dome flats and then a super sky fiat. 

The Keck R band imaging (P.L: G. Illingworth) was 
obtained on November 2003 on a night with photometric 
conditions on the Keck I telescope. 20 exposures of 
300-500 s (average 437.5 s) were ta ken using the re d 
camera on the LRIS Spectrograph ()Oke et al.l [1995D . 
for a total exposure time of 8750 s. LRIS then had a 
resolution of 0.213"/pixel and a field of view of 6' x 7.8'. 
The images were reduced using standard techniques. 
Each indiv idual image was astrometrically calibrated to 
the SDSS (lYork et al.ll2000f ). The final images were then 
comb ined using the drizzle software (jFruchter fc HooS 
|2002[) to a final image scale of 0.1"/pixel. 

The near- infrared J and Ks band imaging (P.L: A. 
Gonzalez) was obtained in December 2003 at the KPNO 
2.1 m telesco pe with the FLAMINGOS instrument 
()Elstonl [19981 ). FLAMINGOS has a resolution of 
0.606"/pixel and a field of view of 20' x 20' on the 2.1 m. 
The imaging was reduced using standard IR imaging 
reduction techniques, using the DIMSUM package of 
IRAF scripts. The final J (resp. Kg) band stacked 
image was made from about 200 frames of 120 s each 
(resp. 890 frames of 30 s each) and has a total exposure 
time of about 24.0 ks (resp. 26.7 ks). 

The Spitzer/mkC (jFazio et al.l [l998l) [i.Q^im\ and 
[4.5//m] band imaging of the clusters (resp. groups) was 
obtained in April 2004 (resp. November 2005 and May 
2006) in 30 exposures of 200 s each (resp. 12 exposures 
of 100 s each) for a total exposure time of 6000 s (resp. 
1200 s) (PI: S.A. Stanford). The data were reduced using 
standard Spitzer procedures. The BCD frames were first 
corrected for muxbleed and pulldown using custom IDL 
scripts (now available from the Spitzer Science Center), 
and then processed with MOPEX to produce co-added 
mosaics. 

2.1. Sample selection 

Our Lynx sample is an ETG subsample of the cluster 
and group homogeneo us sample described in Mil and 
iMei et al.l (|200l [20091 ). ETGs ha ve been visually classi- 
fied in the B rest-frame (zsso) from lPostman et all ([2005[ ) 
for the cluster sample and from Mil for the group sample 
using the same criteria, up to zgso = 24 mag, the limit of 
reliable visual morphological c lassifi cation quantified by 
simulations in [Postman et al.1 ([2005[). At this magnitud e 
our ACS sample is complete (e.g.. [Giavalisco "eFal][200l . 



We apply a selection in photometric redshift (0.92 < 
Zphot < 1.36) and magnitude (21 < zsso < 24 mag). 
The Zphnt. are est imated with Le Phare [3 ([Arnouts et al.l 
[20021 : [Ilbert et~l . 2006) and the selection criteria is cali- 
brated on spectroscopic confirmed members (see Mil for 
details). The magnitude cut at zgso = 21 ensures that no 
star is included in our sample and the cut at zgso = 24 
secures a reliable morphological classification. Galaxies 
belonging to the clusters and groups a re then identified 
by a Friend-O f-Friend algorithm (E OF, [Geller fc Huchra[ 
[1981 see also [Postman et al][2005[ ) with a linking scale 
corresponding to a local distance of 0.54 Mpc, normal- 
ized to 2: = 1.26 a nd to our magnitude range as in 
[Postman et all (|2005[ ) . We use the X-ray emission center 
for the cluste r s and the overdensity centers defined by 
[Nakata et all (f2005l ) for the groups. Spectroscopically 
confirmed outliers were excluded from the sample. 

Our CDF-S sample has been selected following sim- 
ilar criteria. We use as a starting catalog the public 
GOODS-MUSIC v2 sample (ISantini et al.l [20091) . which 
is complete at ^gso = 24 (90% complete at zgso = 26) 
and contains about ~15;000 objects. Photometric red- 
shifts are available for all objects and spectroscopic red- 
shifts, collected from public surveys, are available for 
about 2,900 of those objects. We apply to this CDF- 
S sample the same magnitude cut in zgso band and 
we select objects with secure z^pec (quality fiag=0,l) 
within 1.1 < Zgpec < 1-4. By comparing our selection 
with the number of objects with 21 < zgso < 24 and 
1-1 < Zphot < 1-4, we estimate that our CDES sample 
is more than 70% complete. ETGs were identified by 
visual morphological classification in the zgso bandpass, 
consistent with the Lynx classification. We also verify 
that these ETGs are field ETGs, i.e. that they do not 
belong to already identified struct ures in the CDF-S. Ac- 
cording to lSalimbeni et all (|2009D . there are twelve iden- 
tified structures within the redshift range 0.4-2.5, four of 
which lie at z ^ 1.1. We exclude from our CDF-S sam- 
ple one ETG, which most likely belongs to one of those 
structures. We check the consistency of our m orphology 
classification with that of iBundv et al.l ([2005D : the two 
classifications agree on all galaxies that are in common 
(■^850 < 22.5, 6 galaxies). We thus obtain 27 ETGs in 
the CDF-S with {zspec) = 1.241 ± 0.082. 

The Lynx cluster, group and CDFS-S field samples 
have similar spectral coverage and are complete at zgso = 
24 mag, thus providing an homogeneous and consistent 
sample. A possible bias that might affect the CDF-S 
sample would be a lack of low-mass/passive ETGs that 
are not included in the spectroscopic sample because of 
their faint absorption lines. We will discuss this in our 
result section. 

Our final sample consists of 79 ETGs comprising 31 
in the Lynx clusters, 21 in the Lynx groups and 27 in 
the CDF-S. Our CDF-S sample and about half of our 
Lynx sample have spectroscopic redshifts. We remark 
that known AGNs have been removed from the sample. 
In a companion paper, (Rettura et al., in preparation 
- hereafter Rll), we study in detail the star formation 
histories of the subsample of 13 massive [M > 5-10^°A'/q) 
spectroscopically confirmed ETGs in the Lynx clusters. 



http:/ /www. stsci.edu/hst/acs/ 



http; / /www. oamp.fr/pcoplc/arnouts/LE_PHARE. html 
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TABLE 1 

Lynx cluster images 



Image 


Observation 


Telescope 


Instrument 


Exposure 


FWHM 


Completeness " 




date 






time (ks) 


(arcsec) 


(mag) 


R 


Nov. 2003 


Keck 


LRIS 


8.75 


~ 0.7 


26.2 


«775 


Mar.-Apr. 2004 


HST 


ACS 


7.29 


~ 0.1 


26.3 


^850 


Mar.-Apr. 2004 


HST 


ACS 


12.22 


~ 0.1 


25.9 


J 


Dec. 2003 


KPNO 


FLAMINGOS 


24 


~ 1.6 


21.2 


Ks 


Dec. 2003 


KPNO 


FLAMINGOS 


26.7 


~ 1.4 


21.1 


[3.6/^m] 


Apr. 2004 


Spitzer 


IRAC 


6 


~ 1.6 


23.2 


[4.5/^mj 


Apr. 2004 


Spitzer 


IRAC 


6 


~ 1.6 


22.9 



50% point-source completeness for detection at 5(t 



TABLE 2 
Lynx group images 



Image 


Observation 


Telescope 


Instrument 


Exposure 


FWHM 


Completeness ^ 




date 






time (ks) 


(arcsec) 


(mag) 


R 


Nov. 1999 


Palomar 


COSMIC 


16.2 


~ 1.5 


24.4 


»775 


Dec. 2005 - Feb. 2006 


HST 


ACS 


6.3 


~ 0.1 


26.2 


^850 


Dec. 2005 - Feb. 2006 


HST 


ACS 


10.5 


~ 0.1 


25.7 


J 


Dec. 2003 


KPNO 


FLAMINGOS 


24 


~ 1.6 


21.2 


Ks 


Dec. 2003 


KPNO 


FLAMINGOS 


26.7 


~ 1.4 


21.1 


[3.6/^m] 


May - Nov. 2005 


Spitzer 


IRAC 


1.2 


~ 1.6 


22.2 


[4.5/^mj 


May - Nov. 2005 


Spitzer 


IRAC 


1.2 


~ 1.6 


21.5 


" 50% point-source completeness for detection at So- 



TABLE 3 

Magnitude system conversion (A = AB 



Vega) 



Image 


A = AB - Vega 




(mag) 


_R (Keck) 


0.19 


R (Palomar) 


0.22 


«775 


0.39 


^850 


0.52 


J 


0.92 


Ks 


1.89 


[3.6/im] 


2.79 


[4.5/jm] 


3.26 



3. PHOTOMETRY 

In this paper we will use AB magnitudes in all band- 
passes. TabIS] gives the magnitude conversion between 
Vega and AB system. For ACS and IRAC, those val- 
ues come from the website f^PI . For ground-based tele- 
scopes, those value s have been es timated using the Vega 
spectrum given by iKuruczl (llQQSrFI . 

While SExtractor provides excellent source detec- 
tion and generally good photometry, it has been 
found (e.g., lGiavalisco ern]l20qllBlakeslee et al.l [2006l: 
iHaussler et alj 120071: iMei et alj 12001 " that there are 
some systematics in SExtractor's photometry due to 
sky overestimation. We find similar systematics for 
fixed aperture photometry : when setting the inner ra- 
dius of the annulus used for the sky estimation (with 
the BACKPHOTO-TYPE keyword set to LOCAL), SExtrac- 

http:/ /www. stsci.edu/hst/acs/analysis/zeropoints 
■^"^ http;/ /web. ipac.caltech.cdu/staff/gillian/cal. html 
^ ^ ftp: / / ftp. stsci . odu / cdbs / grid/k93models / standards / vega_k93 . fits 



tor uses the extension of the source defined by SExtrac- 
tor segmentation map and multiplies it by 1.5. In the 
HST/ ACS images, because of the small PSF FWHM, 
this inner radius is in general smaller than the aperture 
radius and this often leads to overestimating the sky, be- 
cause there is still non negligible light from the galaxy in 
the sky annulus. 

In order to take into account the large range of 
PSF FWHMs spanned by our datasct, we perform 
matched aperture photometry (see below), with an 
aperture radius of 1.5" a nd an aperture cor rection 
out to 7" radius (see also IRettura et al.l [2006[) . The 
aperture radius of 1.5", close to the maximum extension 
of our PSF FWHMs (IRAC), is a compromise between 
maximizing the fiux of the source and minimizing the 
contamination by other sources and sky. The aperture 
correction radius, 7", is also a compromise between 
those two opposing goals. 

To obtain accurate photometry, we first build a mask 
for neighboring objects. To this end, we use masks 
obtained with the software SExtractor (we beforehand 
corrected the cases when SExtractor attributes multiple 
detections to the studied galaxy). We then perform 
aperture photometry: we estimate and subtract the 
sky, count the flux within the aperture radius and, 
eventually, apply the aperture correction in order to 
take into account the flux outside the aperture radius. 



3.1. Sky estimation 

Sky estimation is a key step in photometry, because it 
can significantly change the magnitude (in some cases, 
magnitudes derived with SExtractor's sky estimation 
can differ up to ~0.2 from our magnitudes). Due to 
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the variety of our dataset, we set up our own sky 
estimator. We first mask the objects with elhpses, using 
SExtractor's structural parameters hnearly increased 
by a factor of 1.5. We make a first estimate of the 
sky by taking the median value, called skyo, in a 
3" thick circular annulus. The inner radius is set as 
the maximum between 1.5" (the aperture radius) and 
1.5 times the extension of the source. We check the 
galaxy growth curve thus obtained: if it decreases, 
meaning that we overestimate the sky (usually because 
of residual light from masked close bright objects), we 
reduce the sky value by 10% of skyo; we repeat this step 
until the growth curve flattens. This method allows us 
to detect and correct the cases when the sky has been 
overestimated, without any assumption about the shape 
of the growth curve. We do not correct the cases where 
the sky has been underestimated because their detection 
requires assumptions about the growth curve and they 
are rare. 

3.2. Aperture correction 

We measure the flux within a 1.5"aperturc radius. Wc 
then make an aperture correction out to 7", i.e. wc as- 
sume that all the flux of the galaxy is enclosed in a 7" ra- 
dius circle. Aperture correction is usually done using the 
PSF growth curve. Since the difference between the ob- 
served ETG growth curve and the PSF growth curve can 
be significant for the HST/ ACS images, we estimate this 
difference using simulations. In order to build a growth 
curve, we simulate 1,000 galaxies. The simulated galax- 
ies have t he following c haracteristics representative of our 
sample: a lSersid ()1968[ ) profile with an index riser, an axis 
ratio b/a and an effective radius re to which we attribute 
a normal distribution (n^er = 4±2, b/a = 0.65 ±0.2 and 
re = 0.25" ± 0.15") and a random position angle p.a.. 
We then convolve the simulations with the PSF (created 
from selected stars, normalized at 7" radius) and add a 
Poissonian distribution. 

In Fig. [H wc show (left panels) the difference between 
the PSF and simulated galaxy growth curves for the 
two extreme (smallest and biggest) PSF FWHM in our 
dataset {ACS/ HST zgso (top panels) and Spitzer /IRAC 
[3.6^m] (bottom panels)). We can see on the one hand 
that, for the HST/ ACS band (FWHM - 0.1"), the sim- 
ulated galaxy growth curve is much closer, at small radii, 
to the real galaxy growth curve than the PSF one. On 
the other hand, for the Spitzer /IRAC band (FWHM ~ 
1.6"), the simulated galaxy growth curve is very similar 
to the PSF growth curve. However, the difference in total 
flux between the simulated galaxy and the PSF growth 
curves is similar in the two bandpasses. 

On the right panels of FiG. [1] we show a Lynx galaxy 
growth curve calculated using our sky estimate (thick 
black solid line) and the SExtractor's sky estimate (light 
green solid line): this illustrates how SExtractor can 
overestimate the sky in the HST/ ACS images. Pixel- 
to-pixcl sky variations in the HST / ACS sky value are 
high; which make the growth curve less regular. 

Aperture corrections, using PSF and simulated ETG 
growth curves, are displayed in TABjD Interestingly, we 
observe two effects. First, the difference between the two 
aperture corrections is relatively independent of the im- 
age PSF FWHM: it is roughly a difference of 0.1-0.15 




123455 123456 
Aperture radius (orcsec) 



Fig. 1. — Left panels: comparison between the PSF (blue 
dash-dotted Hnc) and simulated galaxy (red dashed line) growth 
curves; if the PSF FWHM is small compared to the galaxy size 
(HST / ACS), the PSF and the simulated galaxy growth curves 
differ significantly at small radii, though it is not the case if 
the FWHM is greater than the galaxy size {Spitzer /IRAC). 
Right panels: example of a Lynx galaxy growth curve using our 
sky estimate (thick black solid line) and using SExtractor's sky 
estimate (light green solid line); we observe that SExtractor can 
clearly overestimate the sky in the HST / ACS bandpasses. 



magnitude. Second, the uncertainty on the aperture cor- 
rection using the simulated galaxies is equal to 0.08 mag 
independent of the band. In order to understand these 
two points, we study the dependence of the aperture cor- 
rection on Te, riser and b/a. We observe that the aperture 
correction is strongly correlated with re, weakly corre- 
lated with riser and independent of b/a. In order to alle- 
viate the dependence on n^er , we repeat our simulations 
with the Sersic index fixed at riser = 4 Ide VaucouleursI 
()1948D profile. We obtain values similar to the ones in 
the right column of TabIH but with a slightly lower dis- 
persion (0.06 instead of 0.08). 

Wc show in FiGj2]how the aperture correction for sim- 
ulated galaxies - convolved with the PSF - depends on 
the effective radius rg. The left panel shows the values for 
four representative bandpasses of our dataset {HST/ ACS 
zgso: magenta stars. Keck R: green dots; KPNO Ks'. 
red squares; Spitzer /IRAC [3.6//m]: blue triangles). The 
right panel shows the same, but the simulated galaxies 
have a fixed Sersic index {riser = 4). The large empty 
symbols at Te = represent the aperture correction de- 
rived from the PSF growth curve. 

This figure illustrates that the dispersion on the aper- 
ture correction mainly results from the dispersion on re 
for the simulated galaxies and explains why this disper- 
sion is independent on the considered band. Moreover, 
it allows us to understand why the difference between 
the two aperture correction methods (PSF and simu- 
lated galaxies) is relatively independent of the considered 
band: the properties of each image are included in the 
PSF, hence the variation with re of the amount of galaxy 
light outward of 1.5" is independent of the PSF. 

In this work, we choose to use aperture corrections de- 
rived from the PSF growth curve (left column of Tab (4]) 
in order to facilitate the comparison with other studies. 
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Fig. 2. — Aperture eorrection as a function of the effective radius re for simulated ETGs. The simulated ETGs have a Sersic index 
following a normal distribution {rtser = 4 ± 2) on the left panel and a fixed Sersic index (riser = 4) on the right panel. For four 
bandpasses (HST / ACS zgso: magenta stars, Keck R: green disks, KPNO Kg: red squares, IPIAC [3.6/^m]: blue triangles) and for each 
simulated ETG, wo represent the aperture correction - i.e. the magnitude correction corresponding to the flux outward of 1.5"- as a func- 
tion of Te- At re = 0, we show (large empty symbols) the aperture correction corresponding to the PSF growth curve (left column of Tab|4]|. 



TABLE 4 

PSF AND SIMULATED GALAXY APERTURE CORRECTIONS 



Image 


PSF 


Simulated galaxies 




(mag) 


(mag) 


R [Keck) 


-0.07 ± 0.02 


-0.18 ± 0.08 


R (Palomar) 


-0.33 ± 0.05 


-0.46 ± 0.08 


1775 


-0.05 ± 0.01 


-0.15 ± 0.08 


^850 


-0.03 ± 0.01 


-0.12 ± 0.08 


J 


-0.33 ± 0.03 


-0.47 ± 0.08 


Ks 


-0.33 ± 0.03 


-0.47 ± 0.08 


[3.6/^m] 


-0.54 ± 0.02 


-0.67 ± 0.08 


t4.5/im] 


-0.56 ± 0.03 


-0.69 ± 0.08 



Note. — The aperture radius is 1.5" and growth curves have 
been normalized at 7". Throughout all this article, wc use the PSF 
correction. 



The photometric and fitted parameters catalogs in Ap- 
pendix are built with the PSF growth curve aperture 
correction. In §5, we analyze the impact of the choice 
of aperture correction on the fitted parameters (stellar 
mass, age). The aperture corrections derived from the 
simulated galaxies growth curve are given in TabIH 

3.3. Magnitude errors 

In order to estimate the error made on the measured 
magnitudes, we used Monte Carlo simulations. For each 
band and each magnitude in the magnitude range of our 
Lynx galaxies with a step of 0.2 magnitude, we simulate 
1,000 galaxies with the characteristics described in §3.21 
We add simulated galaxies in a real image at a random 
position with no object within a 1" radius circle and 
measure their magnitudes with the same method as de- 
scribed above for the observations. So, for a given band 
and a given input magnitude, we now have 1,000 values of 
simulated measured magnitudes. We then estimate the 
magnitude uncertainty as their standard deviation calcu- 
lated with 3(T-clipping iterations. To this estimated un- 
certainty, we add in quadrature a zeropoint uncertainty 
(0.01 mag for HST, 0.03 mag for KPNO, 0.04 mag for 



Keck and Palomar, 0.05 mag for Spitzer). For the PSF 
aperture correction, we also add quadratically the uncer- 
tainty on the aperture correction (left column of Tab|4]). 
If the simulated galaxies growth curve is used, this uncer- 
tainty is already included in the magnitude uncertainty 
derived above. 

4. SED FITTING 

Our catalog permits us to build SED, from which we 
can derive basic galaxy properties. We fit the galaxy 
SED using stellar population synthesis models to obtain 
stellar masses and stellar population ages. 

4.1. Models used 

We compare the results usi ng composite stellar popu- 
lation sy nthesis mod e ls fro m iBruzual fc Chariot (j2003[ ) 
(BC03), IMarastonI ([200l (MA05) a nd C hariot & 
Bruzual (CB07), with a Salpeter IMF (jSalpet er 1955), 
solar mctallicity, no dust and exponentially declining 
star-formation history tjj (characteristic time SFH t: 
ip{t) oc e^*/'"). We thus have a 3D {SFH r, age, mass} 
grid of synthetic spectra to compare with our measure- 
ments. The fitting parameters and their range are listed 
in TabISI 

Because of the age-metallicity degeneracy and to re- 
duce the number of free parameters, we keep the mctal- 
licity fixed during the fit. In H5.21 wc discuss how results 
change when changing the IMF of introducing dust. The 
choice of a SFH as simple as an exponentially declining 
SFH is justified by the lack of information on rest-frame 
UV emitted light (see hereafter), which prevents to con- 
strain more complex SFHs. 

We note that the stellar mass is the mass locked into 
stars, including stellar remnants (column 7 of *.4-Color 
files for BC03/CB07 models and "M A* total" for MA05 
models). The best-fit age output by the models is the 
time elapsed since the onset of star formation. A more 
meaningful age is the star-formation weighted age, which 
represents the age of the bulk of the stars. If we note 
{t)sFW the star- formation weighted age and t the best-fit 
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TABLE 5 

Parameters used for fitting 



Parameter Settings 

IMF Salpeter - fixed 

MctaUicity Solar - fixed 

Dust No dust - fixed 

SFH r (Gyrs) [0.1, 0.25, 0.4, 0.5, 0.75, 1, 1.5, 2, 3, 4, 5] 

Age (Gyrs) [O.l, age of the Universe at the considered redshift] 

Mass(M0) 10* < M <10^^ 



UBVRi775Zs5oJsKs[3.6fim][4:.5fim\. For this set of 
galaxies, we fit the SED both using the complete set of 
bandpasses and only i?i775Z85o Jsi^s[3.6/im][4.5/im]. We 
then compare the ages and stellar masses we obtain with 
those two fits. As it can be seen in FiG. |3l the ages 
derived adding the UBV bandpasses are mostly consis- 
tent with those derived using the bandpasses sampled in 
this paper, though the uncertainty remains large. We 
expected this result, because the 177,5 and zgso filters 
bracket the 4000 A break and thus constrain the age. 
Moreover, as expected, we observe that the masses are 
not significantly affected by the addition of the UBV 
bandpasses. For both ages and masses, adding the U BV 
bandpasses however reduces the uncertainties in the de- 
rived values. 

5. STUDY OF THE SYSTEMATICS IN THE SED 
FITTING 

5.1. Uncertainties due to stellar population synthesis 
models (BC03/MA05/CB07) 

The SED fitting results depend on the stellar pop- 
ulation synthesis model that was used. Until re- 
cently, the most commonly used models were the BC03. 
iMaraston et al.l ()2006D have shown that the treatment of 
the TP-AGB phase of stellar evolution is a source of ma- 
jor discrepancy in the determination of the stellar age 
and mass of high-z galaxies. The modeling of this phase 
is challenging: its rapid evolution makes it difficult to 
obtain the observations needed to constrain the mod- 
els and the physics. Moreover the TP-AGB phase has 
an important contribution at the redshift of our sam- 
ple. BC03 models seem to underestimate the impact 
of the TP-AGB phase, which explains why the inferred 
ages/masses are higher than those inferred with MA05. 
The underestimation of the TP-AGB phase in the BC03 
prescription makes a modeled galaxy less bright and less 
red in the rest-frame NIR. To fit the observation, we then 
need an older and more massive model. The effect of 
the TP-AGB phase is dominant during the first Gyrs of 
the galaxy: taking it into account thus make the galaxy 
brighter and redder at young ages. As our ETGs are ob- 
served when the Universe is about 5 Gyrs old, this effect 
is important. To better take into account this phase. 
Chariot & Bruzual implemented in their new models 
often referred to as CB07) the results of iMarigo et al.l 
2008D on the TP-AGB evolutionary phase. These CB07 
models are still in a preliminary version, but the inferred 
ages and masses are lower than t hose with BC03 models, 
similar to MA05 (see for instance lBruzualll2007b[) . In or- 
der to understand how each model changes our inferred 
masses and ages and to derive conclusions that arc inde- 
pendent of the models, we fit our data with BC03, MA05 
and CB07 set to the same parameters (see Tab. [5]). 

As an illustration, we show in FiG. |31 for one ETG 
of our sample, an EST /ACS zsso stamp, its SED and 
the best-fit spectrum and the output parameters derived 
with each of the three models. Though the three fits are 
of similar quality, the derived age and stellar mass are 
different. 

We compare in FiG. [5] the masses derived with the 
models. Lynx cluster, Lynx group and CDF-S ETGs 
are in red dots, blue triangles and green stars respec- 
tively. To better understand those figures, we also re- 



age output by the models, we have for an exponentially 
declining star-formation history iJj: 

f^{t-t')-^{t')-dt' t-r + r-e-i 

{t}sFw = t = — ; — — (1) 

J^ijit')-dt' l-e X 

Throughout this work, we will consider the star- 
formation weighted age. 

4.2. Fitting method, uncertainties estimations 

We derive the {SFH r, age, mass} parameters by 
choosing the combination which minimizes the de- 
fined by: 

" / \ ^ 

2 \ ^ / rUi data nT-ijmodel \ /n\ 

^=^1 — — r 

z— 1 ^ ^ 

where the index i denotes the band, 771^ ^ the magnitude 
in this band and ai the associated uncertainty. For es- 
timating the la confidence in t ervals , we use the method 
suggested by iPapovich et al.l ()2001[ ). For each galaxy, 
we simulate 250 sets of SEDs by perturbing the original 
photometry randomly within their uncertainties. We fit 
each set with the above method, thus having 250 values 
of minimum Xmin ■ We take the value Xconf ^'^^ which we 
have Xmin < xlonf foi' 68% of the simulated cases. The 
la confidence interval for the original fit is the portion 
of the surface where x^ < xtonf- 

To test how the lack of the measurement in one band- 
pass affects our final results, we have performed a simple 
test on the subsample of galaxies detected in all band- 
passes. When performing the SED fitting on a subsample 
of magnitudes available for each galaxy, we obtain stable 
results only when at least the ACS and one NIR (J, Kg) 
or IRAC bandpasses are available. From hereafter, we 
only consider galaxies that have photometry at least in 
these bandpasses. We thus exclude 4 galaxies for the 
clusters and 1 for the groups: those galaxies are very 
close to another object, which cannot be deblended in 
any of the NIR and IRAC bandpasses. 

We do not consider o ur SFH r estimates a ccurate, be- 
cause, as mentioned in iRettura et al.l ()2010[) , rest- frame 
UV photometry {Xrest < 200 nm) is needed for precisely 
determining the SFH r. In this article, we will restrict 
our analysis to the age and stellar mass; star formation 
histories for cluster ETGs will be studied in our compan- 
ion paper Rll. 

In order to verify that our range in wavelength is 
suitable to derive dependable ages, we used a test 
sample. We consider the subset of 10 ETGs in our 
CDF-S set, for which the photometry is available in 
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Fig. 3. — Reliability of the estimated ages and stellar masses: for a subset of 10 ETGs from our CDF-S set for which the photometry is 
available in the 10 bandpasses from the U band to the [4.5^m] band, we fit the SED once with using all the bandpasses and once without 
using UBV bandpasses. We observe that the ages and masses derived without UBV bandpasses are overall in agreement with those 
derived with UBV bandpasses within the errors. The error bars show the mean error on the estimated age/mass. 
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Fig. 4. — Example of a SED fitting: in the lower right 
panel is a stamp {HST/ ACS 2850 image) of a Lynx ETG 
and, in each other panel, is plotted in red the observed SED 
{Ri77T}Z85oJ Ks[3.6fim][A.5fj,m]) of this Lynx ETG. For each model 
(BC03/MA05/CB07), the best- fit spectrum and the best-fit 
parameters are reported. 



port (squares with a black o utline) data from the litera- 
ture for galax i es at z ~ 1-2 (iBruzuall (l2007aD: inaKeiita : 
iCimatti et aD (|2008| ): light blue: iMuzzin at al.l (|2009f ): 
yellow). The masses are normalized to a Salpeter IMF 
and ages are conv erted to s tar-form ation w eighted ages. 

The sample of IBruzuall (l2007aD is from iDaddi et ahl 
(|2005l ) and has been studied bv lMaraston et al.l (|2006[ ). 
It includes 7 ETGs (1.4 < Zspec < 2.7) selected t o 
be passive with the BzK criteria (jPaddi et al.ll200"^ . 
The SED fitting {V IRAC) for BC03 and MA05 
models is made with various SFHs, a free mctallicity 
and no dust. For the CB07 models, the SED fitting 
uses a simple stellar population with a solar metallicity 



and no dust. The sample of ICimatti et al.l (|2008f ) is 

composed of 13 (1.4 < Zspec ^ 2, mainly early- type) 
galaxies selected in flux ([4.5/zm]) and passively evolving. 
The SED fitting {B — >• Xrest < 2.5/xm) is made with 
exponentially declining SFHs, solar mctallicity and dust 
(the authors note that the inclusion of dust has a weak 
i mpact on the derived ages and masses). The sample 
of|M 

uzzin et aLl (|2009D includes 34 galaxies selected in 
flux (K) at z ~ 2.3 (about half of the sample has a 
spectroscopic redshift). The SED fitting {U IRAC 
+ NIR spectrum) is made with exponentially declining 
SFHs, solar metallicity and dust. 

Though those data cannot be straightforwardly com- 
pared with ours, because the sample selection and the 
SED fitting procedure are different, their addition in the 
figures illustrates the trend seen with our data. When 
comparing BC03 with MA05/CB07, we observe that the 
mass ratio decreases with age for ages lower than ~1- 
1.5 Gyrs and then increases with age. This trend of 
mass ratio with age reflects the activity of the TP-AGB 
phase, which peaks at '-^1 Gyr. In fact, the TP-AGB 
phase activity is important during the first Gyrs (0.2 < 
t/Gyrs < 2 for a sim ple stellar popul ation) with a peak 
around 1 Gyr fe.g.. lMarastonl l2005^ . It is this peak 
shape of the TP-AGB phase activity that we observed 
on Fig. El when comparing BC03 with MA05/CB07, 
but spread on a longer timcscale due to the extended 
star formation history. As we approach the peak of the 
TP-AGB, MA05/CB07 model galaxies arc redder and 
brighter at parity of mass with BC03. This means that 
brighter galaxies are fitted by models with lower mass, 
up to half of the mass given by BC03. We notice that 
the minimum of the mass ratio M(MA05)/M(BC03) is 
reached at ^1 Gyr whereas the one of the mass ra- 
tio M(CB07)/M(BC03) is reached slightly later at - 
1.3 Gyrs. When comparing MA05 masses with CB07 
masses, we also see a trend with the age (the mass ra- 
tio M(CB07)/M(MA05) decreases when age increases), 
however the explanation of this trend is less straightfor- 
ward. 
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Concerning the age (FiG. O, we observe sim- 
ilar but noisier trends, because of larger typi- 
cal uncertainties. The Age(MA05)/Age(BC03) and 
Age(CB07)/Age(BC03) ratios decrease with increasing 
age until ^1 Gyr and then increase with age; the 
Age(CB07)/Age(MA05) ratio decreases with increasing 
age. Our data do not allow us to determine the ages with 
precision (our typical uncertainty on a ge is ^1-1 .5 Gyrs) . 

While pre vious works (e-S-^ iBruzuall l2007al : 
iCimatti et 311120081 : iMuzzin et al.ll2009( ) already pointed 
out that estimated ages and masses depend strongly on 
the chosen stellar population model, our sample allows 
for a better insight into its dependence on the age, which 
was less clear in the three previous studies because the 
sample was too small or too homogeneous in age. 

The photometric catalog and the fitted parameters cat- 
alog are given in the appendix. 

5.2. Uncertainties due to fixed fit parameters 

We now estimate the uncertainty on the derived ages 
and masses due to the parameters we have kept fixed 
(aperture correction, IMF, dust). The mean effect of 
each parameter are summarized in TABjS) We take as 
the reference photometry that which has an aperture 
correction using the PSF growth curve and a fit with 
a Salpeter IMF and no dust. We report the mean and 
standard deviation of the ratio Age/Age^e/ for the ages 
and of log(M/M,.e/) for the masses. 

Choosing an aperture correction with a a mean simu- 
lated galaxies growth curve instead of the PSF growth 
curve leads to a constant shift in mass of about 0.08 dex 
and no significant change in age. This result is expected: 
as the magnitude difference between the two methods is 
roughly independent of the band, using the PSF growth 
curve for aperture correction reduces the fiux in the same 
manner regardless of the band, thus leaving the colors 
unchanged. 

Changing th e Salpeter IMF to a lChabried (|2003( ) (resp. 
iKroupal ([200l for the MA05 models) one results in a 
constant shift in mass of about —0.25 dex (resp. —0.21 
dex ) and no significant change in age. 

We also checked the infiuence of dust extinction by 
adding dust as a free parameter using < E{B — V) < 
0.4, following the iCardeUi et all (| 19891 ) law. We obtain 
on average E{B - 1/) ~ 0.05 - 0.1, except for the CDF-S 
ETGs fitted with BC03 models where E{B - y) - 0.3. 
Adding dust is another way for BC03 models to match 
red colors of young galaxies due to TP-A GB phase. We 
thus c onclude, in agreement with e.g., iRettura et al.l 
(|2006( ). ICimatti et al. (200^, that the assumption of a 
dust-free model for our ETGs is reasonable. 

6. RESULTS 

In this section, we present an analysis of the ages and 
masses estimated with the three different stellar popu- 
lation models (BC03/MA05/CB07). For each figure of 
FiGlTlfTTl we present three figures corresponding to the 
three models, thus allowing a direct comparison of the 
dependence of the results on the models. 

6.1. Color- color diagram 

In FiGl7]we compare observed galaxy colors to model 
predictions. We plot (iyrs — [3.6/Ltm]) against (^775 — -Zsso), 



where age and SFH separate better for our sample. Lynx 
cluster ETGs are in red dots. Lynx group ETGs are in 
blue triangles and CDF-S ETGs are in green stars. The 
models are represented by lines of different colors for dif- 
ferent SFH r (0.1 < T (Gyrs) < 1) and for ages between 
0.5 and 5 Gyrs. 

Lynx cluster ETGs are less scattered than Lynx group 
and CDF-S ETGs and occupy a region corresponding to 
greater ages and shorter SFH t for all models. Lynx 
group and CDF-S ETGs show ages between 1 Gyr and 
5 Gyrs for BC03 and 0.5 Gyr and 5 Gyrs for MA05 and 
BC07. Again this is due to the different modeling of 
the TP-AGB phase. For the range in (1775 — [3.6/im]) 
color between approximately 2.5 and 3.5, observations 
are reproduced by BC03 with ages on average between 
0.5 and 1 Gyrs older than MA05/CB07. We find that 
MA05 and CB07 predict colors that are closer to the 
observations. 

Given the large uncertainties in our SFH r estimates, 
for an accurate analysis see our companion paper Rll. 

6.2. Stellar mass and NIR rest-frame light 

It is well-known that rest-frame the NIR m agnitude 
correlates with galaxy stella r mass (e.g., iGa^azzi 199^ 
iKauffmann fc Charlotl[l99l . In FigIS] we show our esti- 
mated mass versus the [3.6^m] magnitude. The symbols 
are the same as in FigIT] and the solid line is the best 
linear fit given at the bottom of the plots. 

At z 1.26, the [3.6/im] band probes rest-frame H 
where the emission from old stars peaks. We observe 
that the dispersion in this relation is smaller using BC03 
models. MA05 and CB07 show larger dispersion because 
of the different modeling of the NIR emission from the 
TP-AGB phase. To quantify this dispersion, we mea- 
sure the Icr dispersion of the log{Mfit/M) distribution, 
where M is the mass of the galaxy and Mfn the mass 
corresponding to the linear fit. We find a Icr dispersion 
of 0.10 for BC03 models, 0.22 for MA05 models and 0.13 
for CB07 models. For high-redshift galaxies, this effect 
should be taken into account when using NIR rest-frame 
magnitude as a proxy fo r the galaxy stellar mass (see also 
Ivan der Wei et al.ll2006[ ). 

6.3. Ages and masses 

We now discuss the estimated ages and masses. We 
plot in Fig J9] the normalized distribution of the forma- 
tion epochs for Lynx cluster ETGs (red tilted lines). 
Lynx group ETGs (blue tiUed lines) and CDF-S ETGs 
(green horizontal lines) and in FiGdHlthe normalized dis- 
tribution of the masses. For FigJH we consider formation 
epoch and not age, because CDF-S ETGs are observed 
over a period spanning 0.85 Gyrs (1.1 ^ z < 1.4). The 
formation epoch is obtained by subtracting the derived 
age to the age of the Universe at the observation redshift 
5 Gyrs). 

Concerning the formation epochs (Fig19]), we find sig- 
nificant discrepancies between the model predictions. 
The most significant one is for the CDF-S ETGs for 
MA05 and CB07 models on the one hand and BC03 
models on the other hand: MA05 and CB07 models give 
similar formation epochs whereas BC03 models estimate 
a formation epoch on average ^1 Gyr earlier. Because of 
the lack of modeling of the TP-AGB phase, BC03 models 
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Fig. 5. — The dependence of the estimated mass on the stellar population models (BC03/MA05/CB07). Top: the mass estimated with 
one model against the mass estimated with another. Bottom: the mass ratios as a function of age. Lynx cluster ETGs are the red dots, 
Lynx group ETGs are the blue triangles and CDF-S ETGs are the green stars. For our sample, typical mass uncertainty is 40%- 60% and 
typical age uncert ainty is 1-1.5 Gyrs. W e also report (squares with a black outline) data from literature for galaxies at ^ ~ 1-2 IIBruzuall 
Il2007al) : magenta, ICimatti et"all i2003) : light blue, Muzzin et al. (20^): yellow; masses are normalized to a Salpeter IMF and ages are 
turned into star-formation weighted ages). The M(MA05)/M(BC03) and M(CB07)/M(BC03) mass ratios decrease with age for ages lower 
than ~1-1.5 Gyrs and then increase with age due to the effect of the TP-AGB phase modeling. The mass ratio M(CB07)/M(MA05) 
decreases when age increases, however the explanation of this trend is less straightforward. 



TABLE 6 

Variations of age and mass due to the choice of the different fit parameters (see text). 





BC03 


MA05 


CB07 


Parameter 


Age Mass 


Age Mass 


Age Mass 




(dex) 


(dex) 


(dex) 


Aperture correction 


1.06 ± 0.15 0.06 ± 0.04 


1.11 ± 0.30 0.07 ± 0.07 


1.06 ± 0.19 0.06 ± 0.07 


IMF'' 


1.01 ± 0.10 -0.25 ± 0.03 


0.98 ± 0.13 -0.21 ± 0.05 


1.04 ± 0.15 -0.24 ± 0.05 


Dust 


0.57 ± 0.30 -0.11 ± 0.11 


0.82 ± 0.35 -0.07 ± 0.13 


0.75 ± 0.25 -0.03 ± 0.10 



^ For BC03 & CB07 models: Chabrier IMF; for MA05 models: Kroupa IMF. 



artificially increase the age and the mass of the ETGs. 
Ages and stellar masses derived with MA05 and CB07 
seem to be more reliable. 

We note that the mean formation epoch estimated with 
MA05 models for Lynx cluster ETGs is about 0.7 Gyr 
earlier than the one estimated with CB07 models, while 
the mean formation epoch for Lynx group and CDF-S 
ETGs arc similar when using MA05 or CB07 models. 

To quantify the difference in distribution for MA05 
(resp. CB07) formation epochs in clusters and groups, we 
perform a Kolmogorov-Smirnov test and obtain that the 
hypothesis that the cluster and group formation epoch 
distributions are drawn from the same distribution can 
be rejected at a 2.5a significance level: there is a hint 



that the formation epoch distributions are different in 
the clusters and in the groups for MA05 and CB07 mod- 
els. 

Concerning the masses (FiG llOp . we first observe that, 
for the three models, the distribution and the mean value 
of the masses are similar for Lynx cluster and group 
ETGs. Also, there are significantly less massive ETGs 
in the field than in the groups or the clusters (this is 
not the case for BC03 models that artificially increase 
masses for CDF-S ETGs). This is most likely due to the 
low probability of finding massive E TGs in the field and 
is consiste nt with ot her works (e.g.. Fontana ct al.l l2004 
[2006.: .Scodeggio et al...2009: .BolzoneUa et al. 2010,) . 

In the case of the mass distributions, a Kolmogorov- 
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Fig. 6. — The dependence of the estimated age on the stellar population models (BC03/MA05/CB07). Top: the age estimated with 
one model against the mass estimated with another. Bottom: the age ratios as a function of age. Symbols are the as in Fig[5] For our 
sample, typical age uncertainty is 1-1.5 Gyrs. We observe similar trends as in FlG[5]due to the effect of the TP-AGB phase. This trend is 
though more noisy, because of larger typical uncertainties. 




i„5 - [3.6Mm] (AB) 

Fig. 7. — Color-color diagram: Lynx cluster. Lynx group and CDF-S ETGs are in red dots, blue triangles and green stars 
The models are represented by lines of different colors for different SFH r (0.1 < r (Gyrs) < 1) and for ages between 0.5 and 
error bars in the upper right corner shows typical color uncertainties. 



respectively. 
5 Gyrs. The 



Smirnov test docs not permit to reject the hypothesis 
that the two distributions are different. The probabihty 
that they are drawn from the same distribution is of 89% 
(MA05) and 46% (CB07). The masses are then most 
hkely drawn from the same distribution. 

When using MA05 or CB07 models, the average for- 
mation epoch in clusters is consistently earlier than that 
in the groups and in the field. In the field, this might 
be explained by the lack of massive ETGs: galaxies in 
the field are in average younger because they also are 
in average less massive. However, in the groups galaxies 



are in average younger than in the clusters even if their 
mass distribution are most likely drawn from the same 
distribution. 

In FigJTTI we plot the age as a function of the mass, 
with the same symbols as in FicJTl We mark with a black 
outlin e CDF-S ETGs with emission lines (Santin i et al.l 
l2009f ). We observe that, regardless of the model and of 
the environment, age increases with mass, which is con- 
sistent with the downsizing scenario (ICowie et al.l [r9961. 
According to this scenario, the most massive galaxies 
have formed first at high redshifts. This scenario has 
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long seemed to be in disagreement with the hierarchical 
scenario in which m assive galaxies assemb le their mass 
gradually. Recently, iDe Lucia et al.l ()2006D have shown 
that these two scenarios are in agreement if we distin- 
guish the formation epoch of the stars and the assembly 
epoch of the galaxy. Old stellar populations in massive 
galaxies may be assembled by merging of less massive 
galaxies that were already dominated by old populations. 
If however, the CDF-S sample lacks of low- mass/passive 
galaxies (see discussion of the sample selection in H2.1\i 
because of our selection (as it can be seen in FigJTI] the 
low- mass CDF-S ETGs mostly present emission lines ), 
we might miss this population in our analysis. When 



using MA05 and CB07 models, Lynx cluster and group 
ETGs cover the entire range while the CDF-S sample 
shows a lack of old (> 3 Gyrs) or massive (> IO^^Mq) 
ETGs (see above). 

7. DISCUSSION AND CONCLUSION 

In this work, we have studied a sample of 79 red ETGs 
spanning three different environments (cluster, group 
and field) at z ~ 1.3, combining observations of one of 
the most distant superclusters, the Lynx supercluster, 
and of the GOODS/CDF-S field. We built a photomet- 
ric catalog using 1.5" radius aperture photometry with 
PSF corrections, after exploring the possibility of using 
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Fig. 10. — Normalized distribution of the masses: Lynx cluster ETGs (upper, red tilted lines), Lynx group ETGs (middle, blue tilted 
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Fig. 11. — Age as a function o f the mass: the symbols are the same as in FigEI We mark with a black outline CDF-S ETGs with 
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a growth curve which is closer to the real ETG growth 
curve. 

We then estimated the galaxy ages and stellar masses 
through SED fitting using different stellar population 
models (BC03, MA05 and CB07). We show that the 
mass ratio between the masses estimated with differ- 
ent models depends on age. When comparing MA05 
and CB07 with BC03, we observe that their mass ratios 
M(MA05)/M(BC03) and M(CB07)/M(BC03) decrease 
with increasing age until an age of ~1 Gyr and then 
increase. This shape is due to inadequate modeling of 
the TP-AGB stellar phase activity in BC03. Due to this 
problem, BC03 models artificially increase the age and 



the mass to fit the redder and more luminous emission of 
the TP-AGB phase. When comparing CB07 with MA05, 
we also see a variation of the mass ratio with the age (the 
ratio M(CB07)/M(MA05) decreases with increasing age) 
but its explanation is less clear. Concerning the age ratio, 
we sec similar trends with age, though the dispersion is 
greater. The already published mass and age estimations 
on similar ETGs at similar redshifts (IMaraston et al" 
2006 1 lBruzualll2007al: [Cimatti et all 120081: iMuzzin et al. 
2009[ ) with those three models are in agreement with this 
analysis. The advantage of our sample when compared 
to previous results is the larger range in age, that per- 
mits us to identify parameter variations as a function of 
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age. 

This means that current stellar population models give 
uncertain predictions. In the next years it will be es- 
sential to compare their predictions to local and high 
redshift observations to understand their limitations in 
our interpretation of the data. When interpreting ob- 
servations at z ~ 1.3, while some results are stable and 
independent of the model, others depend significantly on 
the different modeling of the TP-AGB phase. 

This inability of the current stellar population 
models to consistently interpret the observations 
has been discussed in various previous works (e.g., 
Maraston et al.ll2006l:rCimatti et al.ll2008l: lEminian et alj 



20081 iLonghetti fc Saraccd 120091: iChen et al.l 
Conrov fc GunnI 120101: IKelson fc HoldenI 
Kriek et al.l " 12011 . IConrov fc GunnI (I2010D 



2010; 



2010; 



have 

compared BC03, MA05 and their own stellar population 
model to local observations and concluded that they 
cannot reproduce star cluster colors in the nearby 
Universe nor the properties of re d sequence massive 
galaxies. Both those authors and iKriek et al.l (|2010f ) 
found that BC03 models better reproduce observed 
SEDs for post-starburst galaxies in local and high 
redshift observations, respectively. 

Our sample at z 1.3 shows that ETG colors and 
SED fits arc slightly better with the new models by M05 
and CB07 when comp ared to BC03 (e.g.. lMaraston et"al] 
I2006t lEminian et al.l [20081 ) . with BC03 predi cted colors 
being bluer than the observations (see also iMei et aTl 
[2006. .2009, ). All models though are consistent with color- 
color diagrams within the uncertainties and give reason- 
able SED fits. What changes is mainly the parameters 
given by these fits and the observation interpretation, 
consistent with what is expected from the new imple- 
mentations of the TP-AGB phase. The advantage of our 
sample is the greater coverage in galaxy ages that per- 
mits us to identify the dependence of the mass and age 
ratios on age, and consequently on the weight of this 
new implementation on these parameters as a function 
of estimated galaxy age. 

Keeping in mind the dependences on the models, we 
obtain the following results: 

1. Independent of the stellar population model and 
the environment, the most massive ETGs show 
consistently older ages. 



According to the MA05 and CB07 models, the 
mass distribution is similar in the clusters and in 
the groups, whereas in our field sample there is 
a deficit of massive (M > IO^Mq) ETGs. This 
lack of massive ETGs in the fiel d , firmly estab- 
lished in the local Univ erse (e.g.. iBell et al.l [20031 : 
iKauffmann et al.ll20Q4[ ). is in agreement at high- 
redshift with recent s tudies in larger samples (e.g. , 
Fontana et all [200l [2001 : iScodeggio etall [20091: 



of massive ETGs. This small difference is less sig- 
nificant than the age difference for galaxies of dif- 
ferent mass (see FiGllip. These res ults are consis- 
tent with previous works at z ~ 1 (jGlemens et al.l 
20091: ICooper et al.l I2010D . Recent results from 



Moresco et al.l (|2010l) (see also lThomas et al.ll20ia 



for similar results in the local Universe) also show 
a stronger dependence of galaxy age on mass and a 
much less sign i ficant dependence on environment. 
IMoresco et al.l (j20100 for example find a difference 
in age of < 0.2 Gyrs, that is consistent with our re- 
sults considering the uncertainties on age estimates 
and that their sample probes less dense environ- 
ments than ours (e.g. docs not include massive 
clusters at z 1). 

In a ACDM cosmological model, galaxies have formed 
first in high density regions and assembled over time 
along filaments to build up larger structures such as 
galaxy groups and then clusters. This cluster assem- 
bly is predicted to happen mainly between 1 < z < 2. 
Our results at z ~ 1.3 show that cluster and group 
galaxies have already formed the bulk of their massive 
ETG population and have a similar mass distribution, 
while in the field we do not observe massive ETGs. This 
might be due to the fact that in a ACDM, massive ha- 
los form in high density regions and/or with a different 
stellar formation history and stellar mass assembly (e.g., 
iPoggianti et~al[ 120061 |2010[) . Regardless of the model 
and of the environment, the more massive galaxies are 
also the oldest, which is co nsistent with the downsizing 
scenario (jGowie et al.lll996[ ). in which the most massive 
galaxies have formed first. This scenario can be recon- 
ciled with the ACDM hierarchical scenario in which mas- 
sive galaxies asse mble their mass from the merging of less 
massive galaxies (jDe Lucia et al.|[2006() if the stellar pop- 
ulations form first in less massive progenitors and then 
assemble later to form massive galaxies. 
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3. When using MA05 and CB07 models, although the 
mass distribution is similar in clusters and groups, 
on average, ages in groups are slightly younger 
0.5 Gyrs). Our field sample population is on 
average slightly younger than that in the clusters 
and groups, consistent with and because of the lack 
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TABLE 7 

Lynx cluster ETG astrometry and magnitudes 



ID 


R.A. 

(J2000) 


DEC. 
(J2000) 


R 

(AB mag) 


«775 Z850 J 
(AB mag) (AB mag) (AB mag) 


Ks 
(AB mag) 


[3. 6^lm] 
(AB mag) 


[4.5fim] 
(AB mag) 


Lynx Cluster E {{z) = 1.261) 



4945 


08 


48 


49, 


.99 


+44 


52 


01, 


,78 


24, 


.81 ± 


0. 


11 


23. 


,68 


± 


0, 


,09 
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,56 


± 


0, 


.07 


21 


.21 ± 0.13 
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19.77 ± 0.08 


19.87 ± 


0.08 
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08 


48 


55, 


.90 


+44 


51 


54, 


,99 
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.58 ± 


0, 


17 
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,30 


± 
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,12 
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,30 


± 


0, 
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U 
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0.09 
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08 


48 
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51 
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0. 


14 


24, 


,20 


± 
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± 


0, 
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U 
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0.09 
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53 
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± 


0, 
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21 
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0.10 
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50 
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33.04 


+44 


53 


39.75 


25.38 




0.14 


24.19 


± 


0.11 


23.30 


± 


0.09 


U 




21.23 


± 


0.14 


20.72 ± 


0.08 


20.92 ± 


0.09 


2094 


08 


48 


34.08 


+44 


53 


32.32 


24.86 


± 


0.11 


24.00 


± 


0.10 


23.13 


± 


0.08 


U 




20.94 


± 


0.13 


20.23 ± 


0.08 


20.24 ± 


0.08 


2343 


08 


48 


35.98 


+44 


53 


36.12 


24.04 


± 


0.09 


22.90 


± 


0.07 


21.89 


± 


0.07 


20.64 ± 


0.10 


19.75 


± 


0.08 


19.17 ± 


0.08 


19.31 ± 


0.08 


2195 


08 


48 


36.17 


+44 


54 


17.30 


24.08 


± 


0.09 


23.03 


± 


0.07 


22.15 


± 


0.07 


20.98 ± 


0.12 


20.08 




0.09 


19.46 ± 


0.08 


19.63 ± 


0.08 


2571 


08 


48 


37.08 


+44 


53 


34.05 


25.07 


± 


0.13 


23.82 


± 


0.09 


22.90 


± 


0.08 


21.52 ± 


0.15 


20.74 


± 


0.10 


20.29 ± 


0.08 


20.48 ± 


0.09 



Note. — Magnitudes are 1.5" radius aperture magnitudes; aperture corrections are estimated with PSF growth curves normalized at 7" radius (cf. 
Tab. magnitude errors are estimated via simulations (cf. section |3.3D ; all magnitudes are in the AB system and are corrected for Galactic extinction. 
U: the galaxy is not detected in the image; B: the galaxy is blended with a close object; -: the galaxy is not on our images. If one wants to use 
aperture correction derived from simulated galaxies growth curve, one has to use the following relations: magGAL — niagpsi^ - ApCpsF + ApCcAi 
and CTQAL — 'sJ^'^sF ~ ^\pC psf^ where magpsi;- and (Tpsf denote the magnitudes and their uncertainty using the PSF growth curve for aperture 
correction (in this table), magG^^, and ctgal the magnitudes and their uncertainty using the simulated galaxies growth curve for aperture correction, 
ApC* the aperture correction and (TApC.psf the uncertainty on the PSF aperture correction (all given TabRI . 
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TABLE 8 

Lynx group ETG astrometry and magnitudes 

ID R.A. DEC. R 2775 2850 J Ks [3.6iim] [4.5tim] 

(J2000) (J2000) (AB mag) (AB mag) (AB mag) (AB mag) (AB mag) (AB mag) (AB mag) 



Lynx Group 1 ((z) = 1.266) 



518 

1339 

1024 

825 

1249 

1085 



08 49 03.52 
08 49 08.32 
08 49 09.00 
08 49 11.24 
08 49 12.27 
08 49 13.69 



+44 53 21.62 
+44 53 48.32 
+44 52 44.08 
+44 51 29.19 
+44 52 13.05 
+44 51 18.82 



24.49 ± 0.14 
23.92 ± 0.12 
24.55 ± 0.15 
23.41 ± 0.10 
B 



23.81 
22.94 
23.72 
22.57 
23.40 
23.04 



0.09 
0.07 
0.09 
0.07 
0.08 
0.07 



23.07 ± 0.08 
21.88 ± 0.07 
22.70 ± 0.07 
21.55 ± 0.07 
22.41 ± 0.07 
22.09 ± 0.07 



U 

20.81 ± 0.11 

U 

20.70 ± 0.10 
B 

21.00 ± 0.12 



22.10 ± 
19.68 ± 
21.22 ± 
19.78 ± 
B 

20.17 ± 



0.33 
0.08 
0.14 
0.08 

0.09 



21.31 ± 0.08 
19.31 ± 0.08 
21.00 ± 0.08 
19.52 ± 0.08 
B 



22.70 ± 0.15 
19.31 ± 0.08 
21.02 ± 0.09 
19.51 ± 0.08 
B 



Lynx Group 2 ((2) = 1.262) 



1636 

939 

1383 

2000 

2519 

1791 



08 49 00.92 
08 49 03.43 
08 49 03.99 
08 49 07.15 
08 49 08.66 
08 49 10.25 



+44 58 49.15 
+44 56 38.59 
+44 57 23.37 
+44 57 52.04 
+44 58 43.26 
+44 56 34.50 



24.45 ± 0.14 
23.64 ± 0.11 
24.70 ± 0.15 
23.78 ± 0.11 
24.86 ± 0.18 
23.17 ± 0.10 



23.45 ± 0.08 
22.87 ± 0.07 
23.65 ± 0.09 
22.91 ± 0.07 
24.09 ± 0.11 
22.33 ± 0.07 



22.35 ± 0.07 
21.80 ± 0.07 
22.68 ± 0.07 
21.93 ± 0.07 
23.18 ± 0.08 
21.42 ± 0.07 



21.28 
21.01 
22.16 
20.71 



U 

20.72 ± 



0.13 
0.12 
0.27 
0.10 

0.10 



20.41 ± 0.10 
20.49 ± 0.10 
20.95 ± 0.13 
19.83 ± 0.08 
21.38 ± 0.15 
19.93 ± 0.09 



19.98 ± 0.08 
20.06 ± 0.08 
20.13 ± 0.08 
19.22 ± 0.08 
20.83 ± 0.08 
19.66 ± 0.08 



20.12 ± 0.08 
20.28 ± 0.08 
20.31 ± 0.08 
19.36 ± 0.08 
20.84 ± 0.09 
19.94 ± 0.08 



Lynx Group 3 ((2) = 1.264) 



137 08 48 53.26 +44 44 22.39 



542 

1135 

889 

1431 

1064 

1136 

1775 

1731 



23.82 ± 0.09 



08 48 55.14 
08 48 56.28 
08 48 56.63 
08 48 57.31 
08 48 57.79 
08 48 57.96 
08 49 01.62 
08 49 04.43 



^44 44 58.83 
^44 46 45.62 
H44 45 39.90 
^44 47 08.01 
^44 45 57.51 
^44 46 04.53 
^44 46 28.23 
^44 45 08.65 



23.26 
23.13 
24.58 
23.86 
23.96 
23.51 
23.52 
22.70 



0.08 
0.08 
0.14 
0.10 
0.10 
0.08 
0.08 
0.07 



22.80 ± 0.08 
22.21 ± 0.07 

22.17 ± 0.07 
23.57 ± 0.10 
23.00 ± 0.08 

23.18 ± 0.08 
22.52 ± 0.07 
22.71 ± 0.07 
21.88 ± 0.07 



U 

20.97 ± 
21.10 ± 

U 

21.77 ± 
U 

21.45 ± 
21.54 ± 
21.19 ± 



0.12 
0.12 

0.18 

0.15 
0.15 
0.13 



21.03 ± 0.14 
20.29 ± 0.10 
20.21 ± 0.09 
21.47 ± 0.18 
21.05 ± 0.14 
21.86 ± 0.26 
20.35 ± 0.10 
20.77 ± 0.10 
20.33 ± 0.10 



20.45 
19.79 
19.77 
20.95 
20.52 
21.42 



± 
± 
± 
± 
± 
± 

19.97 ± 0.08 
20.24 ± 0.08 
19.96 ± 0.08 



0.08 
0.08 
0.08 
0.08 
0.08 
0.08 



19.89 ± 0.08 

19.91 ± 0.08 

19.92 ± 0.08 
B 

20.71 ± 0.09 
21.58 ± 0.10 
20.07 ± 0.08 
20.44 ± 0.09 
20.19 ± 0.08 



Note. — Cf. Tab[7] 
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TABLE 9 

Lynx cluster ETG stellar population ages and stellar masses 



BC03 MA05 CB07 

ID Age log ^ SFH T Age log -M- SFH r Age log ^ SFH r 
(Gyrs) (Gyrs) (Gyrs) (Gyrs) (Gyrs) (Gyrs) 



Lynx Cluster E ((z) = 1.261) 



4945 
6229 
6090 
5355 
8713 
5817 
5634 
5693 
5680 
5794 
8495 
5748 
5689 
5876 
5602 
8662 
8041 
8625 
7653 
8047 
7475 



4.7^ 



7+0.3 
-1.5 

2.7+2-0 



4.0 



0.6 
-1-0.3 
2.2 



1 q-l-1-3 



2.1 
3.0 
4.0 



-1-2.0 
0.8 
+ 1.5 
-0.5 
+0.3 
2.2 



4.4i?:i 

-,+2.0 
-0.3 
a+0.3 
-1.3 
,+0.8 
-0.2 
a+0.3 
-1.0 
1 2+1-2 
1 9+0.6 

a+1-5 
-1.3 
-+0.8 
-0.3 
,+2.2 
-0.5 



1.9: 

s.sj 

2.1^ 
3.8j 



2.8Z 
1.5j 
2.3"* 



11.481°;;^ 
10.90lO;J« 

10.68j:0;^i 
10.65±0;i« 

11 11+0-21 
-'--'-•-'--'--0.05 

11.40+0-?2 



11.44 
10.68 
11.42 
10.98 
10.96 
10.76 
10.43 
10 
10.86 
11.09 



+0.05 
0.18 
+0.16 
0.09 
+0.07 
0.09 
+0.12 
0.07 
+0.03 
0.11 
+0.26 
0.04 
+0.12 
0.05 
+0.13 
-0.12 
+0.09 
0.15 
+0.22 
0.15 



1+0.4 
-0.0 
HO. 7 
-0.0 
s+0.3 
-0.7 
f-0.4 
-0.0 
i+0.8 
-0.3 
HO. 5 
-0.1 

0.810;? 



o.r 

0.1^ 
0.8j 
0.1^ 
0.8j 
0.3^ 



0.1 
0.1 
0.8 
0.1 
1.0 
0.1 
0.3 
0.5 
0.3 
0.3 



+0.4 
0.0 
+0.7 
-0.0 
+0.3 
0.4 
+0.3 
-0.0 
+0.5 
0.3 
+0.3 
0.0 
0.1 
0.1 
+ 0.5 
0.4 
+0.3 
0.0 
+0.8 
0.1 



4.7 
1.9 
2.1 
3.3 
2 

4.3 
1.9 



+0.2 
0.8 
+2.8 
0.5 
+2.6 
0.7 
+0.8 
2.8 
+0.8 
-3.1 
+0.2 
3.1 
+2.8 
0.5 



4 7+<3-2 
' -1.3 
o O+0.8 
0.0_2 g 

o O+0.2 
■^■«-3.8 

0.8t0-6 
0.6+0-^ 



1 

0.6 
1.1 



0.1 
+2.8 
-1.0 
+0.1 
0.1 
+3.6 
0.2 



,+0.04 
-0.07 
-,+0.33 
^-0.16 

10.8llO;f7 

7+0.08 
'-0.43 
10 fiS+0-09 

11.1910-0^ 
11.0810-2^ 



11.433 
10.72"* 



10.77"' 



-, , 42+0-03 

11-*^_0.15 
„+0.09 
" -0.45 
,+0.03 
-0.41 
-+0.14 
-0.36 
,+0.03 
-0.56 
-,+0.17 
-0.15 
,+0.11 
-0.08 

10.59iO-3J 
in 47+0-02 

10 72+0-46 

±U. gg 



10.76: 
11.41^ 
11.05j 
10.93j 
10.60j 
10.13j 



0.1 
0.1 
0.1 


1.5 
0.5 
0.1 



+0.4 
0.0 
+0.7 
0.0 
+0.7 
0.0 
+0.3 
-0.7 
+0.5 
1.3 
+0.3 
0.4 
+0.7 
0.0 



0.1 
0.8 
0.8 
0.1 
1.0 
0.1 
0.1 
0.3 
0.1 
0.1 



+0.4 
0.0 
+0.3 
0.7 
+0.3 
0.7 
+0.7 
0.0 
+0.5 
0.9 
+0.1 
0.0 
+0.1 
0.0 
+0.8 
0.1 
+0.1 
0.0 
+0.9 
0.0 



1+0.5 
1.3 

.+1.6 
'-0.4 

1+1.5 
-0.1 
1+3.5 
-0.3 
!+2.2 
1.2 
>+2.5 
*-0.7 



3.9Z 
1.1 



0.9! 

o.gj 

2.3j 

1., 



3.3 
0.9 
3.5 
1 

3 

0.9 
0.7; 
1.2 
0.8 
1.5 



+1.0 
1.3 
+2.5 
0.1 
+0.3 
2.2 
+1.0 
-0.6 
+0.3 
-3.1 
+0.3 
0.1 
+0.2 
0.1 

+3.3 
0.2 
+0.1 
0.1 
+ 1.3 
0.4 



11.37! 
10.62"* 



7+0.11 

-0.16 
-,+0.25 
-0.23 

10.83tO-i| 
10.35tO-2| 
in 23+o-'"5 

10.94^0-22 

11 ni+0-33 

ll.Ul_Q 24 



n 26+0-1** 

10.33^0-3^ 

-+0.06 
-0.29 
-,+0.18 
-0.24 
-+0.08 
-0.46 
10 64+0-03 

10 23+0-10 
^+0.52 
-0.06 

10.60t0-08 

10.8lt0-26 



11.35! 
IO.8OI 
10.85"* 



10.34"' 



0.4 
0.1 
0.4 
0.1 
0.4 
0.3 
0.3 



+0.1 
0.3 
+0.4 
0.0 
+0.4 
0.3 
+0.4 
0.0 
+1.6 
0.1 
+0.5 
0.1 
+0.5 
0.1 



0.5 
0.1 
0.8 
0.3 
1.0 
0.1 
0.1 
0.3 
0.1 
0.3 



+0.3 
0.4 
+0.7 
0.0 
+0.3 
0.5 
+0.3 
0.1 
+0.5 
0.8 
+0.1 
0.0 
+0.1 
0.0 
+0.8 
0.0 
+0.1 
0.0 
+0.3 
0.1 



Lynx Cluster W ({z) = 1.273) 



1745 


2 


6li 


8 



10 


78^0 


16 
10 





4+0 
^-0 


4 
3 


3 


0+1 
^-2 


3 
5 


10 


761° 


13 
49 





5+0 

O-O 


3 

4 


1 


51? 



3 


10 


43lg 


40 
14 





q+0 


8 
1 


1486 


3 


6to 


3 



10 


7910 


06 
19 


1 


0+0 
^-0 


5 
5 


3 


q+0 


6 

2 


10 


701° 


11 
28 


1 


0+0 
^-0 


5 
6 


2 




8 
7 


10 


5412 


18 
34 





0+0 

^-0 


8 
5 


1794 


3 





6 


11 


03lg 


08 
20 





5+0 

0-0 


3 
4 


2 


3+2 


3 



10 


801° 


25 
21 





q+O 
•^-0 


5 
1 


1 


n + 3 
^-0 


5 
3 


10 


5112 


60 
03 





1+0 


7 



1922 


1 




3 
3 


11 


141° 


12 
02 





1+0 


4 



1 


0+2 
^-0 


5 
8 


11 


031° 


29 
23 





q+0 


5 
1 


1 


5+1 
0-0 


3 
3 


10 


9112 


26 

04 





q+0 
■^-0 


3 

1 


1525 


2 




8 



11 


0410 


17 
10 





4+0 
^-0 


4 
3 


2 


0+2 


6 



10 


85lS 


31 
26 





q+0 


5 
1 


1 


5+2 
0-0 


6 
3 


10 


6912 


38 
15 





q+0 


5 
1 


1962 


3 


Oil 


3 
3 


10 


9710 


12 
12 





5+0 
0-0 


3 
4 


2 


1+2 


6 
7 


10 


79lg 


26 
18 





1+0 


7 



1 


0+2 


2 
3 


10 


6912 


26 
14 





q+0 


5 

1 


2094 


4 


0+? 


3 
8 


11 


211^ 


13 
11 





0+0 

^-0 


3 
7 


4 


1+0 


2 
1 


11 


1712 


14 
37 





0+0 

**-0 


3 
7 


3 


0+0 

^-3 


3 
2 


11 


1612 


06 
59 





0+0 

^-0 


3 
7 


2343 


3 


8+? 


5 
3 


11 


65iS 


11 

13 





c+O 
0-0 


3 
4 


4 


0+0 

^-2 


3 

5 


11 


6210 


10 
31 





5+0 
O-O 


3 

4 


2 


^tl 



8 


11 


4612 


18 
15 





q+0 


5 

1 


2195 


4 


0+? 


3 

5 


11 


54lS 


03 
12 





0+0 

^-0 


3 
4 


2 


ot^ 


6 
7 


11 


211? 


31 
13 





q+0 


5 
1 


2 


411 


8 
5 


11 


3112 


19 
37 





4+0 

*-0 


4 
3 


2571 


2 


s+i 

^-0 


8 

5 


11 


14lS 


17 

07 





0+0 

■^-0 


5 
1 


3 


4li 




8 


11 


1510 


15 
26 





4+0 

^-0 


4 
3 


1 




8 
3 


10 


9112 


32 
11 





1+0 

-'--0 


7 




Note. — Ages {{t)sFw), stellar masses and SFH r arc derived by fitting galaxy SEDs with BC03/MA05/CB07 
models. For the photometry, we use 1.5" radius aperture photometry and aperture corrections are estimated with 
PSF growth curves normalized at 7" radius. If one wants to use the stellar masses derived with photometry using 
simulated galaxies growth curve for aperture correction, see Tab[6] The SFH r are given as a reference for age 
estimation: these values are very uncertain as stated in the text. 



RAICHOOR ET AL. 



TABLE 10 

Lynx group ETG stellar population ages and stellar masses 



BC03 



MA05 



CB07 



ID 



Age 
(Gyrs) 



SFH T Age 
(Gyrs) (Gyrs) 



log 



M 



SFH T Age 
(Gyrs) (Gyrs) 



M 

Me 



SFH r 
(Gyrs) 



Lynx Group 1 ((^) = 1.266) 



518 


1 O+3.0 
1-^-1.0 






^-°-0.4 


0.4 


1339 


■'•■^-i.o 


11 59+°-^^ 




4.0 


1024 


1 9+0.6 
J^-^-0.2 


in fi4+0.l3 


3+°-l 


0.6 


825 


2 1+1-2 
^•1-1.2 






0.8 


1249 










1085 


n 1+2.6 
^•1-0.2 


ll.28t«;,i? 


"■'--0.0 


1.4 




0.7 
2.6 



0.9 
1.8 



+0.4 
0.4 
+ 1.8 
1.7 
+0.1 
-0.1 
+ 1.0 
0.1 

+2.8 



10 16+°-l° 

1U.1D_Q j^g 
3+0.19 
-0.42 
10 ^0+°-"^ 

iU.OU_Q Qg 

™+0.15 



11. 43"* 



11.06" 



0.06 



ll.iO_g 24 



0.3 
0.4 
0.1 
0.1 

0.3 



+0.3 
0.1 
+0.4 
0.3 
+0.1 
0.0 
+0.3 
0.0 

+0.8 



Lynx Group 2 ((z) = 1.262) 



1636 

939 

1383 

2000 

2519 

1791 



9 c: + 2.0 

^•"-•-o.e 

,+0.6 
-0.4 
-,+ 1.3 
-1.3 
-,+0.3 
-1.8 
o O+0.8 
■'•■'-1. 8 
,+0.8 
-0.4 



1-2: 

3.0j 
4.0"* 



1.2^ 



11.24; 
10.99 



+0.16 
0.08 
+0.12 
0.03 
+0.14 



11-18-0.15 

ll.6ll0;f. 
10.97t«;f5 

-+0.14 
" -0.07 



11.15"' 



0.3 
0.3 
0.5 
0.8 
0.8 
0.3 



+0.5 
0.1 
+0.1 
0.1 
0.3 
0.4 
+0.3 
0.7 
0.3 
0.7 
+0.3 
0.1 



9 q + 2.3 
^-■^-1.3 
;+0.1 
-0.1 
1 4+3-3 

,+3.8 
0.1 
+0.1 
0.1 



0.6"* 



0.9 
0.5 



11 11+0-24 
11-11-0.33 

10.7110-03 
10.84lO-« 

11-0'-0.41 

1+0.09 
-0.03 



lO.Sl"" 



Q O + 0.5 


1.8 


1+0-1 
"-1-0.0 


0.7 


0-itro 


1.5 


0.810-3 


3.5 


n 1+0-9 


1.2 


n 1+0-1 

"-1-0.0 


0.7 



s+2.0 
-0.6 
7+0.1 
-0.1 
- + 1.3 
-0.4 
-+0.5 
-3.0 
1-2.8 
-0.2 
1-0.1 
-0.1 



n m+o-23 

il.Ul_g 25 

lu. (y_g 03 
10.791°?^ 
11 'ii+o-io 

11-01-0.51 

in 47+0-41 

1*J-*' -0.03 
10 Q7+0-07 

iu.y<_Q 08 



0.3 
0.1 
0.3 
0.8 
0.3 
0.1 



+0.5 
0.1 
+0.1 
0.0 
+0.3 
0.1 
+0.3 
0.7 
+0.8 
0.0 
+0.1 
0.0 



Lynx Group 3 {{z) = 1.264) 
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Note. — Cf. Tab[9] 
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